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Abstract —Full-duplex small-cell relays with multiple antennas 
constitute a core element of the envisioned 5G network archi¬ 
tecture. In this paper, we use stochastic geometry to analyze 
the performance of wireless networks with full-duplex multi¬ 
antenna small cells, with particular emphasis on the prohahility 
of successful transmission. To achieve this goal, we additionally 
characterize the distribution of the self-interference power of 
the full-duplex nodes. The proposed framework reveals useful 
insights on the heneflts of full-duplex with respect to half-duplex 
in terms of network throughput. 

Index Terms —Full-duplex, multiple antennas, small cells, per¬ 
formance analysis, stochastic geometry. 

1. Introduction 

The employment of full-duplex (FD) technology in wireless 
networks has been recognized as a promising solution to 
cope with the ever-growing demand for high data rates. A 
key drawback of FD wireless communication is the self¬ 
interference received at the FD nodes from their own transmis¬ 
sion. However, thanks to recent advances in self-interference 
mitigation techniques 0^ it is now possible to implement FD 
radios in practical settings In this regard, small-cell (SC) 
systems prove especially suitable for the deployment of FD 
technology due their low transmit powers and the low mobility 
of their users @- 

Several recent works, such as 0, 0, have examined the 
performance of hybrid FD/HD large-scale networks, despite 
considering only single-antenna nodes; in addition, the self¬ 
interference channel gain has been modeled as a constant 
value, which is a very coarse approximation and is only 
meaningful when digital cancellation is applied 0, 0. On 
the one hand, it is timely and relevant to investigate FD nodes 
with multiple antennas in view of the promising concept of 
FD multiple-input multiple-output (MIMO) relays 0, even if 
the extension from the single-antenna to the multiple-antenna 
case sensibly complicates the analysis. On the other hand, 
the residual self-interference channel is known to be subject 
to Ricean fading and, therefore, its modeling in a MIMO 
context represents a challenging problem when receive com¬ 
bining and transmit beamforming techniques are employed. 
In addition, the precise knowledge of the distribution of the 
self-interference power is essential for a rigorous system-level 
performance analysis. 

’Before applying active cancellation, the magnitude of the self-interference 
channel can be modeled as a Ricean distribution with large iT-factor due 
to the strong line-of-sight component; after applying active cancellation, the 
line-of-sight component is reduced, resulting in smaller if-factor 
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Figure 1. System model with full-duplex MIMO small cells ), 

half-duplex backhaul base stations (BHhd), and half-duplex user terminals 
(UThd). with corresponding desired transmissions and interferences. 

In this paper, we fill these gaps by providing the following 
contributions; i) using powerful tools from stochastic geom¬ 
etry, we study the performance of wireless networks with 
randomly distributed FD MIMO nodes and derive tight bounds 
for the probability of successful transmission; ii) we charac¬ 
terize the distribution of the FD self-interference power under 
Ricean fading for arbitrary receive and transmit beamforming 
strategies. In our setting, the FD nodes resemble small-cell 
relays, which are envisioned to be at the foundation of 5G 
©■ Finally, numerical results are reported to corroborate our 
theoretical findings and to establish under which conditions 
the employment of FD is beneficial for the network. 

IT System Model 

A. System Setup 

Consider the scenario in Figure 0 where a set of FD SCs 
acts as relays between a set of backhaul (BH) base stations 
(BSs) and a set of mobile user terminals (UTs), both operating 
in half-duplex (HD) mode; all communications occur in the 
same frequency band. In our setting, the SCs are equipped with 
multiple antennas, whereas the UTs have a single antenna; on 
the other hand, a multiple antenna BH BS performing space 
division multiple access (SDMA) and sending one stream to 
the SC can be considered, hence being equivalently seen as 
single-antenna BS by each SC. During a given time-slot, we 
assume that each SC communicates exactly with one desired 
BH BS and with one desired UT. The analysis of the downlink 











C. SINR Characterization 


(i.e., from the BH BS to the mobile UT) and of the uplink (i.e., 
from the mobile UT to the BH BS) are equivalent since they 
both consist of a succession of a single-input multiple-output 
and a multiple-input single-output transmissions: therefore, 
in the following, we generalize our model and refer to HD 
transmitting and receiving nodes. 

Let us thus introduce the stationary, independently marked 
Poisson point process (PPP) — {(a^n, w(a:„),m(a:„))} 
on X X We use $ = {Xn} to denote the PPP of 
the FD nodes with spatial density A; likewise, $ = fh{^) = 
{fh{x)}x^<s> and <I) = to($) = {m{x)}xe<s> are the isotropic 
marks of $ denoting the HD transmitting and receiving nodes, 
respectively, with hxed distances of the desired links given by 
llcc — fh{x)\\ — R and ||x — fh{x)\\ — R,'ix C <h: evidently, 
$ and $ are PPPs dependent on $ and have also density A. 
In our scenario, it is reasonable to assume that i? ^ ii, since 
the SCs cover a rather small area compared to the range of 
the BH BSs. For convenience, in the rest of the paper we use 
the notation rh^ — th{x) and fhx = fh{x). 


B. Channel Model 


We assume that the FD nodes and the HD transmitting 
nodes transmit with constant power p and p, respectively; 
furthermore, the FD nodes are equipped with Nn receive 
antennas and Nt transmit antennas. The propagation through 
the wireless channel is characterized as the combination of a 
pathloss attenuation and a small-scale fading. The pathloss 
function between the nodes x and z given by l{z,x) = 
\\z — with pathloss exponent a > 2. 

Given transmitting node x and receiving node z, we use 
the following notation]^ The channels are denoted as H^z G 
i^NrxNt if XjZ G $, as hxz G if a: G $ and x G T*, a^ 
hxz G if a: G $ and z G and as hxz G C if a: G 
and 2 G T>; in particular, Hxx models the self-interference at 
a: G $ resulting from its own transmission. Furthermore, Sx 
represents the data symbol transmitted by x with = 

1, whereas the additive noise is denoted by G if 

z G $ and by riz G C if z G $, with elements distributed 
independently as CAf{0,a^). Lastly, G indicates the 
receive combining vector applied by z G $ and w^; G 
is the transmit beamforming vector applied by x G $, with 
llv.f = ||w.f = 1. 

We assume that all the channels, except the self-interference 
channel, are subject to Rayleigh fading with elements dis¬ 
tributed independently as CAf{0, 1). On the other hand, the 
self-interference channel is subject to Ricean fading 0 and, 
therefore, the elements of Hxx are distributed independently 
as CAf{p, In this regard, one can measure the Ricean K- 
factor and the self-interference attenuation H and determine 
the mean and standard deviation of Hxx as (cf. ©) 
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( 1 ) 


In this section, we characterize the signal-to-interference- 
plus-noise ratio (SINR) for the FD nodes and the HD receiving 
nodes, which are needed in the next section to analyze the 
probability of successful transmission. In doing so, we study 
the first hop (i.e., the SINR at the FD nodes) and the second 
hop (i.e., the SINR at the HD receiving nodes) separately under 
the assumption that the whole point process is reshuffled after 


the hrst hop (see details in Section III i. 

First Hop: Consider a reference FD node indexed by k. 
Building on Slivnyak’s theorem m Ch. 8.5], we assume that 
this is located at the origin and, due to the stationarity of $, the 
statistics of its signal reception are seen by any FD node: we 
can thus write l{k,x) = |Xa;|““, with \Xx\ being the distance 
of X from the origin. Hence, the received signal at FD node 
k (its desired transmitter being m^) is given by 
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where (a) represents the desired signal, (b) and (c) indicate 
the interference coming from FD node £ and its associated 
HD transmitting node mg, respectively, and (d) represents the 
self-interference. Given the receive combining vector v^, the 
resulting SINR reads as 


SINRfc ^ ^ 


^rrikk 


( 3 ) 


where we have defined 
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and where Ik is the overall interference at k, i.e.. 


44 {p\x,\- 
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( 4 ) 


^Sek + ^Xff,^\-^Sff,,k)+pSkk. ( 5 ) 


The success probability of the first hop is derived in Sec¬ 
tion UlLA] 

Second Hop: Consider a reference HD receiving node indexed 
by fhk- Again, following Slivnyak’s theorem, we assume that 
this is located at the origin and, due to the stationarity of 
$, the statistics of its signal reception are seen by any HD 
receiving node: we can thus write l(jhk,x) = \Xx\~°‘. Hence, 
the received signal at HD receiving node rhk (its desired 
transmitter being k) is given by 


yfhk= \fpR "h^m^WfcSfc + 

(a) 




+ ^ ^ \/l^Xmt I ^ hjnifhk 


^m£ 




(6) 


2To improve readability, in the rest of the paper, x always refers to a where (a) represents the desired signal, (b) indicates the 
transmitting node and z to a receiving node (either FD or HD). interference coming from FD node I, and (c) is the interference 




















from HD transmitting node m^. The resulting SINR is given 

by 

OTMT3 A ^kfhk 

bINRjJifc — 2 


( 7 ) 


where we have defined 
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^xrrik — 
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and where is the overall interference at rhk, i.e., 

/a, ^ (9) 

The success probability of the second hop is derived in 
Section IIII-BI 


III. Success Probability 

The successful transmission of a packet over the complete 
path, i.e., from the HD transmitting node to the HD receiving 
node through the FD node, is given by the joint ccdf of SINR^ 
and SINRj^j., which is denoted by Pguc — P(SINRfc > 
6>, SINRjfij, > 9 ) pT} for a given SINR threshold 9 (without 
loss of generality, we consider the same SINR threshold for the 
two hops). Assuming no correlation between the two hops due 
to independent sampling of a point process being reshuffled 
after each hop, i.e., the transmission in the two hops occurs 
over two uncorrelated instances of ^rn, it follows that 

Psuc = (10) 

where we have defined P^ic — P(SINRfe > 9 ) and Psuc — 
P(SINRj 7 i^ > 9 ). The case where the interference is spatio- 
temporally correlated will be analyzed in a longer version 
of this paper. Note that, based on the Fortuin-Kasteleyn- 
Ginibre inequality, the uncorrelated case considered here can 
be shown to be a lower bound on the network performance. In 
addition, for notational simplicity, we assume that the system 
performance is interference-limited and, hence, we consider 
4 » 0-2 and Ifh^ > 0-2. 

A. First Flop 

In this section, we analyze the success probability of the hrst 
hop Psic = P(SINRfc > 9 ), i.e., the probability of successful 
transmission from HD transmitting node fhk to FD node k. 
First, considering SINR^ in 0, we have ~ X 2 Nr 

(desired signal) and S^k xi’ G (‘b\{fc}) U ($\{mfc}) 
(interferers)l^ The following lemma gives a tight approxima¬ 
tion of the distribution of the self-interference power Skk- 

Lemma 1. The self-interference power Skk = 

is Gamma distributecj^ with shape parameter a and scale 

parameter b given by 

^ 4i022i4’ ^ 2i2 


^We define a X 2 n r^^idom variable to have pdf f(x) = ^ —. 

^We define a Gamma random variable with shape parameter a and scale 

Of ^ 1 “ X / b 

parameter b to have pdf f(x)= ^ —■ 


respectively, where p and v are the mean and standard 
deviation, respectively, of the self-interference channel H^fc 
(see Q) and where we have defined 

A ‘^NrNt — {Nji + l)(Vr + 1) 

iNji + l){NT + l) 

Proof: See Appendix |I-A| ■ 

Lemma [T] represents a key result of this paper since it pro¬ 
vides a formal characterization of the self-interference power 
experienced by a FD MIMO node with arbitrary beamforming 
vectors, based on the knowledge of the parameters K and H 
(whose values are available either by design or by measure¬ 
ments). 

The next theorem provides the success probability of the 
hrst hop. 


Theorem 1. The success probability of the hrst hop is given 


p(l) A 

^ sue / ^ 
n =0 


z! ds' 




( 13 ) 
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where 

= (i + g5p)a exp(-AT(g)) (14) 

is the Laplace transform of 4 (cf @), where we have dehned 
T(s)^^ (27r-^^^^T'(s,r)ydr (15) 

with 

T'(s,r)= / —(16) 

Jo 1-I-sp(i?2-I--I-2i?r cos 2 

Proof: See Appendix |I-B| ■ 

Given the integral form of T(s) in (m, the success prob¬ 
ability Piic is not in closed-form and needs to be evaluated 
numerically; nonetheless, we derive the following lower and 
upper bounds. 


Corollary 1. The Laplace transform of 4 in ( fT4l l is bounded 
as L/Js) G with 

- (1+7° °’“’ < ^ 

where we have dehned 

T(--)(s) 4 (1 + l){p + p) ^' 7 ^ (19) 

asin(^j 

T(”^-)(s)4 2(p + p) (20) 

asin(-) 

Then, the lower and upper bounds of the success probability 
of the hrst hop Psic are obtained by substituting (s) in 
( |24| with and respectively. 

Proof: See Appendix |I-C| ■ 


Remark 1. In order to efficiently compute the derivatives of 





















the bounds ([T7|i-(fr8|), one can resort to the well-known general 
Leibniz rule | [T2] Eq. 3.3.8] for the differentiation of the prod¬ 
uct of two functions f{s)g{s): for instance, for we 

can write f{s) = and g{s) = exp ( - AT(““)(s)). 

In turn, the derivatives of g{s) can be computed using Faa di 
Bruno’s formula (T^ for the differentiation of the composition 
of two functions {giog 2 ){s), with gi{s) = exp(s) and g 2 {s) = 
_^T(™ax)(g) xjjese considerations apply equivalently to the 
bounds provided in Corollary]^ 

The following corollary provides a sufficient condition 
under which FD outperforms HD in terms of throughput in 
the case of single-antenna nodes. 

CoroUary 2. Assume that = Nt = 1 and let s = 
9p~^R°‘. The lower bound for an achievable throughput in 
the first hop is given by 

= (1 “P (- log3(l + 0) 

( 21 ) 

whereas, if the FD nodes operate in HD mode, the upper bound 
for an achievable throughput is given by 

Th™^(s) = exp A(1 -f £)P ^ log 2 (l + C)- 

\ a sin ^ j y 

( 22 ) 

Then, > T^“^(s) whenever the following condi¬ 

tion holds: 

(1 + sbpr < 2exp - A((l - |)p + 2p) • 

(23) 



Figure 2. Success probability: simulations and theoretical bounds for different 
values of the spatial density A, with Q = —80 dB and 0 = 0 dB. 
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exp ( - AT(“‘^^)(s)) 
exp(-AT(™)( s)) 


(26) 

(27) 


with and defined in ( |T9] l and in ( [20l i, 

respectively. Then, the lower and upper bounds of the success 
probability of the second hop Psuc are obtained by substituting 
(s) in ( |24l l with and respectively. 

Proof: See Appendix |II-B ■ 


B. Second Hop 

In this section, we analyze the success probability of the 
second hop Piuc = P(SINRmj, > 9), i.e., the probability 
of successful transmission from FD node k to HD receiving 
node fhk- First, considering SINR^^, in Q, we have Skfhk ~ 
xInt (desired signal) and ^xrhk ~ xl, Vx € {<^\{k}) U 
(interferers). 

The success probability of the second hop is given next in 
Theorem whereas its bounds are provided in Corollary 


Theorem 2. The success probability of the second hop is 
given by 


N'T — 1 

p(2) A V 

^ sue / ^ 


n—0 


(-s)" d” 
n\ ds’ 


(s) 


(24) 


s=ep-^R° 


where 

^T'(s,i?)exp(-AT(s)) (25) 

is the Laplace transform of (cf. (|^), with T(s) and 
T'(s, r) defined in (|T5]l and inlfTh]), respectively. 


Proof: See Appendix II-A| 


Corollary 3. The Laplace transform of in ( |25l l is bounded 
as £/_ (s) e with 

k 771^ 


IV. Numerical Results 


In this section, we present numerical results to assess our 
theoretical findings and, specifically, to compare the perfor¬ 
mance of FD with respect to HD. In the following, the pathloss 
exponent is a = 4, the distances characterizing the marked 
PPP are set to i? = 5 m and R = 0.5 m, the transmit 
powers are p = 0.5 W and p = 1 W, and the considered 
SINR threshold is 0 = 0 dB; the parameters a and b of 
the self-interference power are computed according to o, 
where p and iz are obtained from Q with Ricean AT-factor 
K = 1 (see Q for an experimental characterization of K) 
and self-interference attenuation = —80 dB. Lastly, the FD 
nodes adopt maximum ratio combining and maximum ratio 
transmission and, therefore, the beamforming vectors are given 
by 




(28) 


Figure plots the success probability Pguc in ( flOl t over 
different values of the density A. Remarkably, the theoretical 
bounds obtained in the previous section are reasonably tight. 
Furthermore, the employment of multiple antennas allows to 
mitigate the effect of the self-interference of the FD nodes, 
thus producing substantial SINR gains (observe that to higher 































Figure 3. Comparison between FD and HD modes in the first hop: minimum 
throughput gain for different values of the self-interference attenuation fl, 
with A = 10“^ and 6 = 0 dB. 



Figure 4. Comparison between FD and HD modes in the first hop: minimum 
throughput gain for different values of the SINR threshold 6, with A = 10“^ 
and Q = —80 dB. 


values of Nn correspond more terms in the summation ([T3|)). 

We now focus our attention on the first hop (subject to self¬ 
interference) in order to analyze the benefits of FD. With this 
objective in mind, we introduce the minimum throughput gain 


rp^(mirL) A 


rTn{min) 

^FD 

rTn(max) ’ 

-^HD 


(29) 


with and defined in ( |2T] l and ( |22] l, respectively: 

this parameter denotes the worst-case gain of FD mode over 
HD mode in terms of throughput, with > 1 indicating 

that FD outperforms the equivalent HD setup. Figure plots 
TG(min) over different values of the self-interference atten¬ 
uation n with A = 10“^, whereas all the other parameters 
are the same as in the previous simulation. In this setting, 
we have > 1 even for moderate values of the 

attenuation and, in the specific, when: H < —52 dB for 
Nr ^ Nt = I, n < -44 dB for Nr = Nr = 2, 
n < -36 dB for Nr = Nr = 4, and H < -28 dB for 
Nr = Nr = 8. Moreover, the minimum throughput gain is 
analyzed in Figure as a function of the SINR threshold 9 
with A = 10“^ and fl = —80 dB: in this respect, it is shown 
that FD achieves improved performance with respect to HD 
for any reasonable value of 9. 


V. Conclusions 

This paper analyzes the performance of wireless networks 
with FD MIMO small cells using stochastic geometry. We 
characterize the distribution of the self-interference power 
of the FD nodes for arbitrary receive/transmit beamforming 
strategies and we derive tight bounds for the success proba¬ 
bility. Our framework highlights the beneficial effect of FD, 
which produces substantial throughput gains over HD for 
realistic values of number of antennas, node density, self¬ 
interference attenuation, and SINR threshold. 


Appendix I 

Success Probability of the First Hop 


A. Proof of Lemma [T] 

Due to space limitations, we only provide here a sketch of 
the proof. For notational simplicity, in the following we omit 
the sub-indices in the beamforming vectors and in the channel 
matrix and write Skk — |v^Hwp, with 

NNt 

|v^Hw|2 = ^ ^ v*Vkh,jhlgWjw} (30) 


where Vi is the z-th element of v = Wj is the j-th 

element of w = and hij is the (z,j)-th element of 

H = Recalling that hij ~ v'^), we have 

E(|/iijp) = + andE(|ft,yj^) = /j,^-|-4p,^i/^-|-2iz^; on the 

other hand, for any normalized receive/transmit beamforming 
vector, 

1 


we have E(|wi 


“ Nr^Nh + I) J — 

Nr(Nt+i) • ®*iilding on the central limit theorem for causal 
functions we approximate a sum of positive random 

variables x = Xi by a Gamma distribution with shape 
parameter a and scale parameter b defined, respectively, as 

(E(a:))^ A Var(a:) 


A 

a = 


Var(a;) 


A 

a = 


E(a:) 


( 31 ) 


In order to compute a and b, we need to derive the second 
and fourth moments of jv^Hwp. Let ai denote the z-th 
singular value of H and A^min = xam{NR, Nr): from 
after some algebraic manipulations we obtain 


Npt N'j' 

E(|vHHwH =e( ggkn/z. 


Wi 


i=i j=i 


2 I 2 

— p -f IZ 


(32) 

(33) 





























( 38 ) 


Ci^{s) =E{exp{-spSkk))^i n exp(-s(p|X£| °‘Sek + p\Xfhi\ ^Smik)) 

^ee<s>\{k} 


n exp(-s(p|X^| Smk + p\Xfht\ (41) 

^ f U\ / 


'lG^\{k} 


Nrr 


and 

E(|vHHw| 4 ) =e( J2 \v, 




+ 2E 5] 


A^min 


|f,;PhiP|CTiP|tT,-P|w,;P|w,-P 


* J = 1 
*7^i 


4iV^iVr 


(-^i? + 1)(-^T + 1) 


other hand, the Laplace transform of If^^, is obtained as in ( |4T] ) 
at the top of the page using the moment-generating function 
of the xi distribution. Again, we resort to ||^ Th. 2] and the 
expression in (|25]l readily follows. ■ 


(34) 


/ + + 2i^^. (35) 


B. Proof of Corollary 

Given the definition of 'I'(s,r) in ([T^, we observe that 
1 1 


4'(s, i?) G 


Since Var(|vHHw|2) = E(|vHHw|4) - (E(|vHHw| 2))^ 
we readily obtain a and b in 0 from ( [3 T] i. ■ 

B. Proof of Theorem [T] 

The success probability of the first hop is given by 


1 + sp(i? — i?) “ l + sp(i? + i?)“ 


(42) 


where we recall that i? ^ i?. Then, the lower and upper 
bounds of £/_ (s) in @-(l^ are a straightforward result 
of combining (|4^ and Corollary ■ 
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where Ik is defined in (|^ and denotes the overall interference 
at k. Since FD node k is equipped with Nn receive antennas, 
the power of its desired signal is distributed as S^^k ^ X2 Nh 
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Appendix II 

Success Probability of the Second Hop 
A. Proof of Theorem 

The success probability of the second hop is given by 

p(2)^ = P> 0 \ (39) 

\ ^fhk J 

= E{Skfh,>Op-^R‘^Ifh,) (40) 


where Iff^^ is dehned in (|^ and denotes the overall interference 


at rhk- Since FD node k is equipped with Nt transmit 
antennas, the power of the desired signal of rhk is distributed 

Skfhk 

expression in 


xInt 


and, similarly as in Appendix 


results from applying 115 Th. 1] 
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